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Abstract 


An  expression  for  the  scattered  wave  as  an  expansion  in 
terms  of  radial  eigenf unctions  is  obtained.  The  total  wave  may  be 
expressed  in  terras  of  radial  eigenf\mctions  directly  but  the  in- 
cident wave  cannot.  Instead,  the  incident  wave  is  first  expressed 
by  a  contoiir  integral.  Ihen  a  change  of  variable  is  introduced  and 
the  resvilting  integral  is  approximated  by  the  Eviler-Maclaurin  sum 
rule.  This  resuD.ts  in  a  series  for  the  incident  wave  in  terms  of 
radial  fxinctions  plus  an  integral  and  correction  terms.  Wten  this 
is  svibtracted  from  the  total  wave  a  finite  series  in  terms  of  radial 
functions  is  obtained,  and  the  integral  and  correction  terms  are 
easily  evaluated. 
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1,  Introduction 

This  paper  is  concerned  -with  obtaining  the  amplitvide  of  plane  waves 
scattered  in  the  forward  direction  by  a  sphere  whose  radius  is  large  compared 
to  the  wavelength.  Both  Dirichlet  (vanishing  function)  and  Neumann  (vanishing 
normal  derivative)  boundary  conditions  will  be  treated. 

The  procedure  used  for  this  problem  is  applicable  to  many  electro- 
magnetic and  quantum  mechanical  problems  dealing  with  waves  scattered  by 
spherically  stratified  regions  with  various  boundary  conditions. 

The  description  for  the  scattered  wave  in  terms  of  partial  waves  has 
practical  limitations  for  short  wavelengths  because  of  poor  convergence.  For 
the  total  wave  this  difficulty  has  been  overcome  by  means  of  an  alternate  re- 
presentation characterized  by  an  expansion  in  terms  of  radial  eigenfunctions 
(see  Friedman,  Franz L  J)  rather  than  angvilar  eigenfunctions  which  appear  in 
the  partial  wave  solution.  The  purpose  of  this  paper  is  to  show  a  method  for 
expressing  also  the  scattered  wave  through  the  radial  representation. 

These  two  alternate  forms  for  t  he  total  wave  exhibit  distinctly  dif- 
ferent convergence  properties.  The  expansion  in  terms  of  angular  eigenfunctions 
is  more  rapidly  convergent  for  plane-wave  scattering  and  large  wavelengths.  On 
the  other  hand,  the  radial  representation  is  more  rapidly  convergent  for  small 
wavelengths  when  the  source  is  close  to  the  sphere  (cf.  Marcuvitz  ^ -J ,  p.  310) . 
Since  the  scattered  wave  alone  is  produced  by  sources  (induced  by  the  incident 
wave)  directly  on  the  sphere,  the  radial  eocpansion  promises  to  be  more  rapidly 
convergent  for  the  scattered  wave  for  short  wavelengths. 

The  incident  wave  cannot  be  subtracted  directly  from  the  total  wave 
in  its  radial  form,  because  the  radial  eigenfunctions  are  constructed  to  vanish 
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on  the  surface  of  the  sphere,  in  keeping  with  the  boundary  conditions.  Therefore, 
no  combination  of  these  eigenfunctions  can  duplicate  the  value  of  the  plane  wave 
on  the  surface. 

Instead,  the  incident  field  is  represented  first  by  its  expansion  in 
terras  of  partial  waves;  this  expansion  can  then  be  expressed  as  a  contour  integral, 
by  applying  a  Watson  transformation.  Through  the  use  of  the  Euler-Maclaurin  sum 
rule,  a  portion  of  the  integral  may  be  represented  by  a  series  of  radial  eigen- 
functions. When  the  result  is  compared  vath  the  sum  for  total  wave,  this  radial 
representation  does  indeed  show  the  expected  convergence  properties.  Only  a 
few  terms  in  the  series  for  the  total  wave  needed  be  retained,  and  the  remaining 
portion  of  the  integral  for  the  plane  wave  is  easily  evalxiated. 

The  final  numerical  values  for  the  first  order  frequency  dependent  cor- 
rection to  the  total  cross  section  for  both  boundary  conditions  agree  essential- 
ly with  those  obtained  by  other  authors  (cf.  Wu  and  Riibinow,  White  "-'J),  The  method 
used  here,  however,  is  not  primarily  dependent  on  computational  methods  and  would 
appear  to  have  more  physical  significance. 


2,  Fonnfulation  of  the  problem 

Consider  a  sphere  of  radius  a  whose  center  lies  at  the  origin  of  the 
set  of  spherical  coordinates  as  shown  in  Tig.  1 


Figure  1 
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Incident  on  the  sphere  and  propagating  in  the  positive  z-direction 
is  a  plane  wave  u.  of  wave  nuniber  k: 

/ o  T \  ikr  cos  6 

(2.1)  ^  u.  «  e         • 

The  reduced  wave  equation  that  descidbes  the  total  field  for  the  corresponding 
acoTistic,  electromagnetic,  or  quantum  mechanical  problem  is  given  by 

(2.2)  V^u  +  k^u  »  Of 

for  this  particular  problem  we  specify  the  impedance  boundary  condition 

■^  -  kzu  =  0    for  r  =  a  , 

(2.3) 

^  -^  iku     as  r  -J^oo  . 

Ihe  solution  for  this  problem  expressed  as  an  expansion  in  terms  of 
Legendre  functions  has  been  previously  obtained  (Friedman,  Franz  ••-'): 

CD  A  (ka)    /^\ 

(2.1,)  u    '  22    j„(kr)i°(2n+l)P„(cos  6)  -  T]  — hy'(kr)i''(n-H/2)P  (cos  Q)   , 

n=0  '^  "         n   B^(ka)    "  ^ 

where  A  =  2j  (ka)  -  z  2j  (ka)  and  B^=>  h:""-^  (ka)  -  z  h^-'-^(ka)  .   Note  that 
li         n  n  n   n  n 

for  r  =  a,  equation  (2,U)  satisfies  the  condition  prescribed  by  (2.3), 

Note  also  that  the  first  sum  in  (2,h)  is  the  expansion  for  the  plane 
wave  in  terms  of  the  P  ,  namely  (cf.  Strattonl-J) 

(2.5)         u^  =  e^^  cos  e  ^  ^  j^(kr)i"(2n+l)P^(cos  6)  . 

Therefore  the  remaining  sum  in  equation  (2,li)  represents  the  scattered  wax'Co 
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However,  equation  (2,ii)  is  not  the  onl7  means  of  expressing  the 
total  field.  An  alternative  expression  for  u  may  also  be  written: 


(2.6)  u 


n 


\   (ka) 
n 


1> 

"n-  2' 


{- 


in 


(v„+  i)n  expj-  :^  (v„+  1) 


n 


sin  V  n 
n 


P^  (cos  ©)h^^^(kr)  . 
n       n 


v=v 


n 


In  this  equation,  the  values  v  are  obtained  from  those  values  of  v  that  satisfy 


(2.7) 


B^  S  h^^^  (ka)  -  z  h|,^\ka) 
n     n  n 


0 


Equation  (2.6)  may  be  arrived  at  by  applying  the  Watson  transformation 
(see  i*!riedman,  Franz'--')  to  (2.I4)  or  by  a  technique  using  separation  of  vari- 
ables and  utilizing  the  properties  of  the  Green's  functions  and  eigenf unctions 
of  r  and  (e)(see  MarcuvitzLJ,  pp.  309  ff.). 

We  note  that  h^   (kr)  satisfy  the  radial  equation  (2,2),  that  they 
n 
satisfy  the  boundary  condition  at  r  =  a  as  indicated  by  (2,7),  and  that  they 

are  outgoing  as  r  becomes  large.  Therefore  h   (kr)  are  the  radial  eigen- 

n 
functions,  and  (2,6)  is  the  above-mentioned  alternate  expansion  of  the  solu- 
tion in  terms  of  the  radial  eigenf unctions.  The  convergence  properties  of  these 
two  representations  (2.6)  and  (2,Ii)  have  already  been  mentioned  in  the 
Introduction, 


3,  Subtraction  of  the  plane  wave 

We  wish  to  obtain  a  representation  for  the  scattered  wave  that  will 
exhibit  the  convergence  properties  of  an  expansion  in  terms  of  radial  modes. 
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Since  the  plane  wave  cannot  be  split  off  directly  fl*om  (2,6),  as 
it  «as  from  (2,U),  it  must  be  subtracted  in  some  other  manner.  To  this  end, 
we  first  represent  the  plane  wave  given  by  (2,5)  by  a  contour  integral.  Ac- 
cordingly, we  may  write  u.  in  the  form 


r    3  (l(r)(2v 

(3.1)      "1  =  ^  J r 


2v+l)  e 


ivn/2 


P^(cos  9) 


xvn 
e    sin  vn 


dv 


where  the  residues  reduce  to  the  series  (2,5)  when  the  contour  is  taken  as 
shown  in  Fig,  2, 


/3 


(v  ♦  l/2)-plane 


7 
Z 


Figure  2 


(3.2) 


The  lower  branch  of  the  contour  may  be  reflected  by  introducing 


-  V  -  1  , 


because  then 
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»  V  +1/2=  -00 

(3.3)   I   2J  (^rV^)"  e-^^/2  P  dv  -   I    2j      (^'^1/?)"  ^i(v'+l)n/2 

J  ,  J  sin  V  n  -V '  -1 

v+1/2  -  0  V  +1/2  -  0 

But  in  general 

(3.1.)  24>,e-i""/2  .   ^^1),  ,^2))  ,-1-/2  ^ 

and,  using  the  identities 

-V -1  V         -V-1  V 

and 

(3.6)  P  ,  (cos  9)  =  P  ,(cos  e)     for  all  v   , 

-V  -1  V 

we  may  vrrite  instead  of  (3«U) 

r^  i\  PS  i(v'+l)n/2  -iv'n/2  ..  /^^  "/^  .  (l)   .  ^^^  ^\ 

(.3«7;  2j     ,     e  '     »    e  '     2J  ,  -  e  h^,     2  cos  v  n     , 

-v-1  V  V 


Now  the  integral  in  (3»3)  becomes 

2J 


J        ^Sl4^«-^^"/^  2J^P^dv  .  J     2J   ,     Iv^l4)n     ^-iv'n/2  ,  ^  ,^. 


V        sin  V  n 
v+1/2 -0  v'+l/2  =  0 


(3.8)  -00 

v'+l/2=0 


-00 

.     r     2h^l)     £2£JVi     e^^'"/2(v'+l/2)n  P  ,  dv'. 
J  V        sin  V  n  v 
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The  original  integral  about  the  positive  real  axis  in  the  v  +  ^  -plane 
now  becomes  an  integral  above  the  real  axis  from  +00  to  -00  plus  an  additional  con- 
tour starting  from  the  origin  to  -co.  Removing  the  primes,  we  obtain  for  u. 


(3.9)  V  ^ 


"'+03 


as  illustrated  in  Fig,  3« 


■r 


(v^l/2)n      -ivn/2  p  ^^         f      ^iS^)  cosjm  ei^«/2(^^i/2)„p  ^v 
"^v    sin  vn  v  J         ^      sin  vn  >■     */   '     y 

v+1/2  -  0 


P 


(v  +  l/2)-plane 


Fig\ire  3 


According  to  (2.7),  there  will  be  values  of  v  -  v^  that  specify  the 

zeros  of  B     (ka)  .     Assxime  that  the  v  's  delineate  a  path  C  in  the  v+  ^  —plane 

n 
as  illustrated  in  Fig.  h',  then  deform  the  contours  for  u.   in  the  manner  shown 

by  the  figure. 
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V  +  -k  plane 


Figure  U 


The  integrals  over  the  large  semi-circle  may  be  shown  to  vanish  as  the 
radius  grows  larger.  Hence  the  integral  representation  for  the  incident  wave 
may  be  written  in  terms  of  the  integral  whose  contour  starts  at  the  origin  and 
is  deformed  over  the  path  delineated  by  the  radial  eigenvalues. 

r 


(3.10) 


u^=  . 


^ 


2h^^\^)   g^  e^^"^  (v-.l/2)n  P^(cos  e)dv  . 
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Now  the  scattered  wave  u  is  given  by  the  total  wave  miniis  the  incident 

s 

wave.  We  write  u  as  the  svm  for  u  given  by  (2,6)  minxis  the  integral  in  (3«10) : 


-3-  iZI 

V 

n 


s 

(v^+  l/2)n 

sin 


A^(ka) 


l-B  (ka) 


'  +  l/2)n     r  .„        -I  f^s. 

2j^^^  sxpL-  f  (v„*  1)J  P„  (cos  6)h^l)(.a-) 
n       *-         -J   n        n 


(3.11) 


r 

i 


2h^l^kr)  g|^  ei"^/2  ^^  ^  ^/2)n  P^(cos  e)  dv  . 
C 

In  the  following  sections  we  shall  insert  in  (3.11)  a  change  of  vari- 
able in  both  the  sum  and  the  integral  and  introduce  the  tangent  approximation 
for  A  and  B  in  the  sum.  We  recognize,  though,  that  the  tangent  approximation 
to  the  Bessel  functions  is  not  valid  for  values  of  v  +  1/2  near  ka.  Therefore, 
this  approximation  does  not  apply  to  the  initial  terras  of  the  sum.  After  the 
first  L  terms  however,  when  v,+  1/2  is  sufficiently  far  from  ka,  the  tangent 
apprciximation  becomes  acceptable  and  may  be  used  in  the  remaining  terras  of  the 
sum.  The  criterion  for  determining  L  depends  on  the  accuracy  desired. 

We  intend  to  show  that  the  tangent  approximation  for  the  siun  effectively 
cancels  the  integral  in  (3 .11)  when  its  contour  conforms  to  the  path  delineated 
by  the  eigenvalues.  To  achieve  this,  we  specify  Vj.+  1/2  as  the  point  at  which 
the  contour  begins  to  conform  with  the  delineated  path  (see  Kig.  ii). 

In  later  sections,  the  initial  terms  of  the  sum  are  evaluated  by  the 
more  accurate  Hankel  approximation.  Also,  the  initial  portion  of  the  integral 

(from  V  +  1/2  equals  zero  to  v^+  1/2)  is  easily  approximated  if  terms  of  order 

-6/3 
(ka)  '   are  neglected. 
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ito  Evaluation  of  the  total  wave 

Asymptotic  forms  for  the  splierical  Hankel  functions  axe  obtained  by  the 
saddle-point  method  applied  to  the  Sommerfeld  integral  representation  for  solu- 
tions to  Bessel's  equations.  The  forms  that  are  appropriate  here  are 


(U.l) 


a.2) 


V 


h(«(ka) 


2e 


-inA 


ka 


ka 


-  1 


-lA 


cos  (i-  +  w)  +  ..•  , 


1 


V  + 


-  1 


ka 


-lA 


-iw  ^ 
e    +  ••• 


with  the  condition  that  the  real  and  imaginary  parts  of  v  +  ^  are  both  greater 
than  zero  and  that  the  function 


(a.3) 


w  -  y  ika)^-   (v  +^)^  *  -  (v  +  |)  cos 


-1  (11^' 
ka 


has  a  real  part  less  than  zero.     Equations  (h.l)  and  (U»2)  are  forms  of  the 
so-called  tangent  approximation.     (cf»  Bremmer'    '). 

The  derivatives  of  (U.l)  and  (U«2)  are  needed  also.     These  are  given  by 


(U.U)  h^^^ 


-  2 


ka 


M^-3 


-lA 


dw 


sin(]-  +  w)  TT--  +  a  negligible  term. 


("•«  4^''-  ^hh^J-^ 


-lA 


^-iw  dw    _^   ^  negligible  term 
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For  brevity,  let  us  introduce  4  and  C  defined  by 


(U.6) 


dw 


^^  =  C  sin  4  and  Z  =  C  cos  ?  • 


Then  B  reduces  to 


(U.7)  B^  = 


(:?'- 


znt"^ 


-  2 


C  e 


-in/h 


la 


M'-  r 


cos($  -  T-  -  w)   • 


The  zeros  of  B  occur  when  the  cosine  term  in  (i^.?)  vanishes,  li.   a  similar 


manner  we  find 


(h.8)  A^  =   2j^-  Z  2J^ 


iC 


M-r^^'-<(^i^'"-i^). 


Now  vre  introduce  a  change  of  variable.  Let 
(1.9)  S  -  J  -  ^  "    "(^  +  |) 


Since  the  real  part  of  w  must  be  less  than  zero,  x  is  restricted  to  positive 
values.  In  particular,  when  x  is  zero  or  an  integer,  A  vanishes.  This  cor- 
responds to  V  *  V  provided  the  tangent  approximation  is  valid.   -The  derivative 
of  B  with  respect  to  v  is  found  from  (l;.?)  to  be 


(i^-10)  |7^=  2 


^^  [(^T-  ] 


-lA 


sin  n(x  ♦  |)  ^  " 


+  a  cosine  term  that  vanishes  at  v   • 


The  ratio  of  (U»8)  and  (ii.lO)  at  v  =  v  becomes 


(U.ll) 


57  ^v 


e 


i(x+l/2)n 


„  dx  .  /  ^  1\_ 
Vn      ^"  37  ^"^"^   2^" 


-  1 


2ni^ 
dv 


-Jv=v 


n 
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The  derivative  -r—  is  found  frcm  (lv«3)  to  be 
ov 


(1.12) 


3w 


-  cos 


Therefore,   (Ii.ll)  reduces  finally  to 


(U.13) 


A 

V 


2  cos 


n 


^^^     ' 


Equation  (2.6)  for  the  total  wave  (or  the  sum  in  (3»11))  is  written 

now  as  two  sums.  The  first  sum  extends  from  v  to  v, •  The  tangent  approxima- 

o  Jj 

tion  is  not  acceptable  for  the  terms  in  this  sirni.     The  second  sxim,   however, 
extends  beyond  v^j  for  this  sura  the  tangent  approximation,  and  therefore 
(U.13),  is  acceptable.     Accordingly,  we  write  ti  as 


v^      r- 


u 


L 

V 


A 
V 

i  l-B 
3v     V 


— 'V 


(V  lA)  ^[-  ^  ("„-  1)] 


sin  V     n 
n 


P     (cos  e)  h^^\kr) 
n  n 


n 


(U.IU) 


00 

-ZZ 

V 


L+1 


2  cos 


-1     V  +  1/2 
ka 


(v^+1/2)   exp[-^  (V^O 


sin  V     n 
n 


P     (cos  e)h^.-'-}kr) 


n 
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5,  Evaliiation  of  the  incident  wave 

The  incident  irave  is  given  by  (3»10).  This  is  the  integral  subtracted 
from  the  series  for  the  total  wave  in  (3«ll)» 

Let  us  write  this  integral  in  two  parts,  denoted  by  L.  and  I_.  The 
first  integral,  L ,  extends  from  the  origin  to  the  L-th  eigenvalue,  and  the 
second  integral,  I„,  extends  from  the  L-th  eigenvalue  outward.  These  paths 
are  illustrated  in  Fig.  5» 


/3 


v-i-  ^  plane 


Figure  5 


First  consider  the  integral  1^.     We  introduce  the  change  of  variable 
defined  in  ih»9)   and  (U.3): 


(5.1) 

Then 
(5.2) 


y^(ka)^-  (v  +|)^  '  -  (v  +  |)  cos"^ 


2  -  n(x  +  #)  +  C 


dw«^dv  =  -n^. 
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^  has  already  been  found  from  (U.12) ,  Substituting  this  value  into  (5«2) 
we  have 

n  dx 


(5.3)  dv 


cos 


■^(^ 


Next,  we  use  an  asymptotic  form  for  cos  vn: 


/  r-  \  \  „  1  -ivn  f.         +2ivn  \     1 

(S.U)     cos  vn  =  -  e     f  1  +  e     I   =  ^ 


-ivn 
e    '<'•••  • 


In  a  later  section,  we  will  demonstrate  that  the  imaginary  part  of  v,  is  like 

Im  Vj^  ^  (-3  Wj^)^/^  (ka)^/^  . 

Then,  for  any  value  for  Im  v  along  the  path  taken  by  Ip,  the  asymptotic  terra 
in  (5«b)  is  smaller  than 


Which  approaches  zero  faster  than  any  negative  power  of  ka  as  ka  approaches  in- 
finity. When  ka  is  as  low  as  unity,  neglecting  the  asymptotic  term  results  in 
an  error  of  less  than  1  part  in  20,000  at  the  L-th  eigenvalue. 

Using  (5.3)  and  the  asymptotic  form  of  (5»U)  in  the  integral ( 3,10) ,  we 
may  write  I_  as 


(5.5)   i2=  -I 


'     /-  ^  -.  /o\        -in(v+l)/2  /^x 

^^•^y^)"  .  •  2—, —    P  (cos  e)h^.^^kr)  dx 


cos 
y'x=L 

If  we  consider  the  change  of  variables  introduced  earlier,  i.e.,  com- 
pare (5.1)  with  (i4.3)  et  seq.,  we  see  that  x  is  a  real  variable  along  the  path 
prescribed  for  Ip  and  that  v  takes  on  the  values  v  whenever  x  is  an  integer© 


-  15  - 


Let  MS  now  evaluate  the  integral  (5»5)  by  the  Euler-Maclaiirin  sum  rule 
(see  Whittaker  and  Robinson  ^-l).      This  series  is  the  following  expression  for 


a  definite  integral: 


{S.6) 


Jx«a 


=a+NAx 


f(x)dx 


■[K- 


•) 


n=N 


f(a)  -  f(a+Ni:^)J     +  212   f(a-»TiAK) 


Zix  +  5 


This  expression  is  equivalent  to  truncating  the  curvB  for  f (x)  as  shown  in  i-'ig,  6c 


f(x) 


Figure  6 


across  successive  values  of  /J^  x  and  summing  the  areas  of  the  trapezoids,     llie 
term  5  in  (5»6)  is  a  fcorrection  term  for  the  error  resulting  from  the  trvulcations 
and  is  given  by 


(5.7) 


00 


(-D^'^Ax^"^*^ 


1^0    "^^         (2m  +  2): 


.2m+l 
f (a+NAx)  -       2^^^     f (a) 
dx 
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The  coefficients  B^  -■  are  called  the  Bernoullian  nunibers.  The  first  few  of 
<:in+l 


these  are 
(5o8) 


1 
^1  Z  ' 


B,  =  ?  ,   B^  =  tt?  ,   B,  =  175-  , 


2^ 

30 


1 
^3  "  H? 


>  ••  • 


For  our  purposes,  a  =  L  and  ^x  =  1,  which  corresponds  to  the  interval 
between  the  successive  eigenvalues  v  •  Also  it  na.j  be  shown  that  the  integrand 
in  (5.5)  as  well  as  its  derivatives  vanish  as  N  approaches  infinity.  For  the 
case  at  hand,  then,  the  Euler-Maclaurin  rule  becomes 


(5.9) 


/° 


x=L 


00 


f  (x)dx  =  i  f  (L)  +  211  f  Cn)  +  6  , 

n-L+1 


where  f  is  the  integrand  of  (5.5)»     Accordingly,   5  becomes 


(5.10) 


00 


(-1) 


ra 


,2m+l 


S)    ^'""l     (2ra+2):       "^^ 


f(L) 


Ip  may  now  be  written  as 


(5.11)  Ig  -  ^  f(L)  + 


n»L+l 


(v^+  1/2 )n 


cos 


-i/V^^^" 


-in(v„+l)/2  .X 

5 2 p^  (cos  0)hj,^^kr) 

sin  V  n     n        n 
n 


+  6, 


Corapsirison  with  equation  (U.lh)  for  the  total  wave  shows  that  the  sum 
appearing  in  (5.11)  agrees  exactly  with  that  in  (U.lli)  for  n«L+l  and  beyond.  Thus 


I-  may  be  written  as 


(5.12)   Ij-^Z 


A 


( 


v^+  l/2)n  .  exp[-^  (v^+  1)] 


sin  V  n 

n 


P^  (cos  ©)h^^\kr) 
n       n 


+  u  +  I  f(L)  +  5  , 
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and  the  incident  wave  may  be  expjreseed  by 


V,   r— 


\=    u  +  £3 


(^.13) 


A 


3v  V 


(v^^  1/2)ti  exp[-  ^   (v^H-  1)] 

sin  V  n 
n 


P^  (cos  e)hj^^^kr) 
n       n 


n 


+  1^+1  f(L)  +  5  . 


($.13)  shows  that  the  difference  between  u  and  u.  i.e.,  the  scattered  wave,  con- 
sists of  a  finite  series,  a  definite  integral,  and  a  correction. 


6.  Evaluation  of  the  scattered  wave  in  the  forward  direction 

In  order  to  evaluate  the  definite  integral  I- ,  we  introduce  certain 
restrictions.  First  we  consider  only  the  far  field,  for  which  r  »  a.  Then 

h   (kr)  maybe  written  asymptotically  as  (cf»  Stratton  L-'J ) 

n 


(6.1) 


V      kr 


We  also  restrict  the  problem  to  forward  scattering,  for  which 


(6.2) 


P,(l) 


for  all  V  • 


Then  the  integral  I. ,  which  is  (3«10)  extended  from  zero  to  v,*  1/2,  can  be 

expressed  as 

(vi+  1/2) 


(6.3) 


h 


1   e^^ 
7    '^ 


tan(v  +  i)Ti  •  (v  +  i)n  dim 


This  integral  is  most  easily  evaluated  along  the  path  shown  in  ^±g»   7» 
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i0, 

/ij  1    1 

/"u'     2 

/ 

» 

/ 

i 

• 

/ 

0 

a 

Figure  7 


Let  us  denote  the  Ira(v  +  "t)  by  Pj   then  the  integral  along  the 


imaginary  axxs  xs 


(6.U) 


as  xs  . 

I         tan(v  +  |)n(v  +  |)n  dvn     «     -  i     j       ta 


tanh  pn  pn  dpn     • 


A  graph  of  tanli  pn  •  pn  appears  in  Tig.  8j  the  shaded  area  shows 
how  it     departs  from  a  straight  line.     This  shaded  area  is  finite  no  matter 
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tonh  ^TT  -  ^ir 


Flgurt  8 


how  far  pn  is  extended.  Let  us  call  this  area  A A,  and  define  it  in  terms  of 
(6.U)  as 


n\ 


(6.^ 


-  i  I   tanh  pn  pn  dpn  =  -  i  j   pn  dpn  +  i  i£SA(p^n)   , 


For  that  portion  of  the  integral  extending  from  lfl_  to  v  +  i  shown 
in  Fig,  7,  the  asymptotic  forms  fcr  the  cosine  (and  sine)  specified  by  (5«U) 
may  be  used,  yielding 


tanh(v  +  |)n  =  i  f  1  +  ofe"^^^   jj. 


Therefore,  the  entire  integral  I.  may  be  written 

ikr 


(6.6) 


T  -    i   ® 


I 


(Vj^+  1/2) 


(v  +  i)n  dvn  +  A  A 


which  reduces  immediately  to 
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(6.7) 


h 


=  -  i 


ikr 
e 


^\^ 


1/2)' 


-7 


We  turn  now  to  evaluating  the  finite  series  given  in  (^•13).  To 
this  end,  we  write  first  the  Hankel  approximation  for  h^  ,  2j  ,  and  their 
derivatives  (cf.  Bremmeri--') 

C6.e)    h'l'dca)  .  ^   f^    [j^/3(-w)  .  J.i/3(-w)] 

(6.9)    2J,  (to)  .  1.  /^  [,-i"/3  .^/^(.H)  .  ei"/3  J.^/3(.K)] 


(6.10)   h("'(to)  -  ^  f^   [.,/3(-w)  -  ..,/3(-w)]  .  (S,) 


(6.11)   2^:  (to)  -    ^Y^    p"/3  .,/3(-«)  -  e-i"/^  J.3/3(-w)]  .  (^)   . 

In  these  equations,  w  has  been  given  before  by  (h«3),  and  p  is  given  by 

./    V  -^  1/2  2 ' 

Note  also  that  (6.8)  and  (6.10)  can  have  zeros  only  for  negative  real  values  of  w. 

Let  us  consider  the  Dirichlet  condition  first,  letting  z-^00.  This 
reduces  the  bracketed  term  of  the  sum  in  (5.13)  to 


(6.12) 


i  l-B 
dV  V 


V 

-J  n 


2o^(ka) 


n 
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where  the  v  are  obtained  from  h.  (ka)  =  0.  (6.8)  shows  that  these  roots 
n  V 

n 


follows  from 
(6.3) 

We  need  as  well 


V3^-V  -  -  '^-1/3^-V  • 


(6.11,) 


n 


+  a  vanishing  term  at  v  . 


(6.13),  (6.1ii)  and  (6.9)  reduce  (6.12)  to 


(6.15) 


where 


(6.16)   D. 


A 


-D 


n 


V 

— '  n 


(V3^-V^^:i/3(-ni)(i),^      Pt 


,   siibject  to  (6.13)  • 


We  call  D  the  Dirichlet  coefficient,  where,  according  to  (U.ll),  D  approaches 
w  as  the  tangent  approxiraation  becomes  acceptable. 

The  Neumann  boundary  condition  for  which  z  »  0  is  treated  in  an 
analogous  manner.  For  this  condition 


(6.17) 


A 


^¥\ 


n 


2j^  (ka) 


V 

— '  n 


where  the  v  are  obtained  from  h\   (ka)  »  0.  (6.10)  shows  that  these  are  the 

roots  of 

(6.18)  '^2/3^  ~V  "      '^-2/3  ^"V  • 
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The  derivative  of  (6.10)   reqtdred  for  (6.17)  is 


(6.19) 


+  a  vanishing  term  at  v  • 


(6.11),  (6.18)  and  (6.19)  reduce  (6.17)  to 


(6.20) 


where 


A 


-  N 


n 


V 

— '  n 


f2/3'-V  -  'UM-\^-  {^\  (^ 


n 


(6.21)   N^  =   r'  r        /      r*         r  '  subject  to  J2/3^"V  '     '^-2/3^~V  * 


V3^-^n^  -  '^-2/3^-V 


We  call  the  N  the  Neumann  coefficient,  which  must  also  approach  ^  as  the  tan- 
gent approximation  becomes  acceptable.  Through  the  use  of  (6.1),  (6.2),  the 
asymptotic  form  for  sin  v  n,  and  either  (6.l5)  or  (6«2C),  the  sum  appearing 

in  (5.13)  reduces  to 

ikr         (v^+  1/2) 
(6.22)  -  2ni  V  Z:  C„  7^ 

'vn 


M~' 


where  C  is  either  D  or  N  as  the  case  may  be. 
n         n    n 

Also,  the  correction  term  6  must  be  evaluated.  This  is  obtained  from 
(5«10),  where  the  fxmction  f  is  the  integrand  of  (5.5).  After  inserting  the 
appropriate  asymptotic  forms  sind  specializing  to  the  forvjard  direction,  f  becomes 

-ikr 


in 


(6.23)  f  =  in  -j^ 

In  (6.23),  equation  (ii«12)  has  also  been  used. 


(v  *   1/2) 
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(6.2h) 


In  both  (6.22)  and  (6.23),  the  term 

y  +  1/2 
dv/dv 


appears.  We  now  obtain  an  approximation  for  this  term.  This  is  done  by  in- 
verting ih'}))   which  yields  (cf.  Bremmer^-J) 


(6, 


.25)    v-H  1/2  =  ka  [l  +  I  (ka)-2/^(-3w)^/^  e^"/^  +  0(ka)"^/^]  . 
dw/av  is  fo\md  by  differentiating  (6.25): 
(6.J6)      ^   -  -  [|-3w)-l/3(to)l/3  ,1V3  ,  0^^)-h/iy    . 


Therefore  (6.2U)  becomes 

(6.27)    i^  =  .(ka)2  [(ka)-2/»(-3w)-^/3  ,iV3  ,  0(ka)-''/3]  . 

The  integral  L,  given  by  (6.7),  may  be  expressed  at  once  through 


(6.25)  as 

(6.28)  Ii=  -i-j^ 


ei^  Oca)!  [^1  ,  (ka)-2/3(.3w^)V3,in/3  ,0(ka)-^/3.  2Aa  (^).6/3"j. 


Substituting  (6.27)  into  the  series  given  in  (6.22)  yields 
(6.29) 


2ni  4^  (lca)2  ei"/3 


kr 

Also  (6o23),  again  by  (6.27)  becomes 

ikr 


^    C^(ka)-2/3(-3w^)-l/3  ^  0(ka)-^/3 


(6.30)  f  =  -  in 


kr 


(kH)2  [(k.)-^/3(-3w)-V3  e^VJ  ,  0(ka)-''/3]  . 


The  correction  term  (see  (5.10))  requires  the  derivative  of  f  with 
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respect  to  x  at  v^+  1/2,  This  is  done  through  ih»9) 


(6.31) 


,2!n+l 


dx 


2m+r 


(-3w) 


-1/3 


,2ra+l 


dx 


where  $  =  0  for  the  Dirichlet  condition,  and  n/2  for  the  Neumann  condition*  The 
asymptotic  form  for  the  scattered  wave  in  the  forward  direction  may  be  written 
now,  using  (5»13)  and  the  preceding  equations: 


-ikr 


u  «  u  -  u. 

S  X 


i  e^  O-ll   1  .  (^)-2/3ei"/3 


(-3wj^)2/3  +  (-3wj^)~^/^.  n 


(6.32) 


-  hnT^    C  (-3w  )"^/^  +  2n5-'  B  ^^  -iiL_ 
^   "    "         'fe'o  ""^l  (2m+2): 


0 

,2m+l 


dx 


2S^ 


[3n(x  .  3A)  -  34]""-/^ 


where  higher-order  term  have  been  omitted. 

We  shall  compute  now  the  values  of  (6»32)  for  the  two  boundary  condi- 
tions. Consider  the  Dirichlet  condition  first.  We  compare  the  roots  found  from 
the  Hankel  approximation  (6.13)  with  those  of  the  tangent  approximation,  namely 
the  zeros  of  cos(n/l4  +  w)  in  (U»l) : 


Tangent  approx. 

Hankel  approx. 

•*w 
o 

2,35619 

2,383l;6 

-*1 

5oU9778 

5,51020 

-*2 

8.63937 

8,614736 

-*3 

11.78096 

11,78685 

-25- 


The  last  values  quoted  above  differ  by  about  1  part  in  2000.  The 
values  of  D  (6.16)  corresponding  to  the  roots  from  the  Hankel  approximation 
are  given  in  the  next  table: 


^n 

V 

o 

.h9h9S 

^1 

.h9Q9S 

^2 

.1*9959 

^3 

.U9980 

These  values  clearly  approach  1/2,  the  last  term  differing  by  about 

2  parts  in  5000.  The  correction  terra  is  found  from  the  first  term  of  the 

series  for  6  to  be 

5  =  -  .00709  • 

Higher  corrections  give  contributions  beyond  the  fifth  decimal  place. 

After  the  remaining  arithmetical  computation  is  performed,  we  obtain 


(6.33) 


Ug(©  =  0) 


(ka)2  ["1+  (ka)-2/3  e^"/3x  .9966q\ 


for  the  Dirichlet  condition. 

We  turn  finally  to  the  Neiunann  condition.  First,  we  colI^5are  the 
roots  determined  by  the  Hankel  approximation  (6.18)  with  those  of  the  tangent 
approximation,  i.e.,  the  zeros  of  sin(n/l4  +  w)  in  (Ij.h)? 
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Tangent  approx. 

Hankel  approx. 

-w 
o 

0.785hO 

0,6851|2 

-^1 

3.92699 

3.90279 

-W2 

7.06858 

7.05I49U 

-^3 

10.21018 

10.20068 

The  last  terms  differ  by  about  1  part  in  lOOOo 


The  values  of  N  (6.21)  obtained  from  the  Hankel  s^proximation  are 


fo\md  to  be 


\ 

w 

0 

.5U978 

*1 

.50295 

*2 

.50093 

"? 

.5ooUii 

which  also  approaches  1/2,  but  in  the  opposite  sense  from  the  D  's.  The  cor- 
rection term  is  found  again  from  the  first  term  of  the  series  for  5  to  be 

5  -  -  .00861 


and  the  remaining  conqsutation  results  in 

ikr 


(6,3U)    u_(e  =  0) 


i  nST-  (ka) 


[=■ 


(ka) 


-2/3  ein/3 


.8651 


for  the  Neuaann  condition. 
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The  total  scattering  cross  section  may  be  found  by  multiplying  the 
imaginary  part  of  the  scattered  amplitude  by  l^n/k  (cf»  B ohm  •--»),  The  final 


results  are 


(6.35)  C     =  2n  a^  Tl  +  (ka)"^/"^  .99668  \  Dirichlet  case  , 

(6.36)  a-     =  2n  a^  n.  -  (ka)'^'-^  .865ii  )  j  Neumann  case 

These  results  are  in  substantial  agreement  with  those  found  recently 
by  Wu  and  RiibinowL'Jj  who  obtain  values  of  .99618  and  -.86140  for  the  two 
values  above. 

It  is  of  interest  to  note  that  if  the  tangent  approximation  is  used 
throughout  (which  simplifies  the  procedure  outlined  in  this  report)  an  error 
of  approximately  3  percent  is  foiind  for  the  Dirichlet  condition,  but  consider- 
ably greater  error,  about  30  percent,  occtirs  for  the  Neumann  condition. 

7.  Conclusion 

The  significant  feature  of  the  method  employed  in  this  problem  is 
the  manner  of  subtracting  the  plane  wave.  It  was  shown  that  the  scattered 
wave  could  be  represented  by  a  finite  series  plus  a  definite  integral  (apart 
fi*ora  the  correction  6).  (5.13)  shows  that  this  can  be  done  without  restrict- 
ing the  general  impedance  boxindary  conditions  or  the  values  of  r  and  e.  Also, 
it  can  be  shown  that  the  subtraction  applies  for  the  source  located  at  a  finite 
position;  one  need  not  restrict  the  method  to  an  incoming  plane  wave. 
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It  is  necessary  to  depart  from  generality  only  to  evaluate  the 
definite  integral.  This  would  not  be  necessary  for  the  sum,  since  appro- 
priate approximations  exist  for  P  as  vrell  as  for  h^  •  Therefore,  ex- 
tension to  the  more  general  case  appears  plausible  and  deserves  further 
study* 

The  distinction  between  the  cases  treated  here  and  more  general 
cases  lies  in  locating  the  singiilar  points,  hence  the  radial  eigenvalues. 
But  the  form  and  accuracy  of  the  solution  would  be  comparable. 
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